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The studies of Fermi-Dirac correlations of AA and AA pairs in hadronic Z decays, Bose-Einstein correlations 
and colour reconnection in W-pairs decays performed by the ALEPH collaboration in e + e~ annihilation at LEP 
are presented. The OPAL analysis of Bose-Einstein correlations in W-pair decays is also discussed. 



1. Introduction 

Bose-Einstein (BE) correlations between iden- 
tical bosons and Fermi-Dirac (FD) correlations 
between identical fermions lead to an enhance- 
ment or a suppression, respectively, of the parti- 
cle pairs produced close to each other in phase 
space. The effect is sensitive to the distribu- 
tion of particle sources in space and time. The 
strength of the correlations can be expressed by 
the two-particle correlation function C(p\,p2) = 
P(p 1 ,p 2 )/Po{pi,P2), where pi and p 2 are the 
four-momenta of the particles, P{pi,P2) is the 
measured differential cross section for the pairs 
and Po(pi,p2) is that of a reference sample, iden- 
tical to the data sample in all aspects except the 
presence of FD or BE correlations. Usually C(Q) 
is measured, where Q 2 = — (pi — P2) 2 ■ 

For W-pairs from e + e~ — > W + W~ at energies 
in the LEP2 range, the distance between W + and 
W~ vertices is less than 0.1 fm, i.e. less than the 
typical hadronic distance scale of 1 fm. Therefore 
the fragmentation of W + and W~ may not be in- 
dependent. Two phenomena may appear: pions 
from different W's may exhibit BE correlations 
and pairs of quarks and antiquarks qiq 4 and q3q 2 
from the decay of different W's can form colour 
strings (colour reconnection). Colour reconnec- 
tion (CR) and BE correlations have opposite ef- 
fects. They may influence the accuracy of the W 
mass measurement at LEP. 

In this talk, three ALEPH analyses are pre- 
sented: FD correlations of AA and AA pairs in 
hadronic Z decays at LEP1 M, colour reconnec- 
tion Q and BE correlations |3[ in W-pairs decays 
at LEP2. The OPAL analysis of BE correlations 



in W-pair decays |J] is also discussed here, due to 
the unavailability of the OPAL speaker. 

2. Fermi-Dirac correlations in (AA, AA) 
system 

The FD correlations in (AA, AA) system were 
studied using 3.9 million hadronic Z decays 
recorded by the ALEPH detector on and around 
the Z peak. A sample of 2133 pairs with Q < 
10 GeV was obtained. 

In the analysis, three reference samples 
were used: A) simulated pairs from JET- 
SET MC without FD correlations, where 
C(Q) = P(Q) dat jP(Q)uc; 

B) pairs obtained by event mixing, where 

C(Q) = [P(Qh,tjP(Q)^ a ]/[P(Q)uc/P(Q)7ch 

C) reweighted sample of mixed pairs, where 

C(Q) = P(QWJP(Q)™Z:^ d - . 

The measured correlation functions are shown 
in Fig. [IJ parametrised with 

C(Q) =Aqi + /3exp(-i? 2 Q 2 )] (1) 

Consisted results are obtained for the three ref- 
erence samples. The correlation function C(Q) 
decreases for Q < 2 GeV; as Q tends to zero, 
it approaches the value of 0.5, as expected for a 
statistical spin mixture ensemble. If this is inter- 
preted as a FD effect and the parametrisation of 
Eq. (|l|) is used, the resulting values for the source 
size R and for the suppression parameter (3 are 



R= 0.11±0.02 stat ±0.01 sys 
P = -0.59 ± 0.09 stat ± 0.04 sys 



fm 



An alternative method to study the (AA, AA) sys- 
tem is to measure the spin composition of the 
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Figure 1. Correlation function for the (AA, AA) 
pairs using reference samples A, B, and C. The 
curves represent fits using the parametrisation 
given in Eq. ([!]). 

system using the angular distribution dN/dy*, 
with y* the cosine of the angle between the two 
protons (antiprotons) in the di-hyperon centre-of- 
mass system. The measured dN/dy* distribution 
has contributions from both S = 1 and S = 
states: 

dN dN 
dy* 



dy* 



■ £ 



5=0 

where e is the fraction of 
tion. By fitting the dN/dy* 



dy* 



S=l 

5 = 1 contribu- 
distribution in dif- 



ferent Q ranges, the dependence e{Q) is ob- 
tained. But e(Q) can also be defined as e(Q) = 
C(Q) s =i/[C(Q) s=0 + C(Q) s =i], with C(Q) s =i 
and C(Q)s = o the contributions of S = 1 and 
S — states to the correlation function C(Q). 
Using the parametrisation given in Eq. ([!]), one 
obtains for a statistical spin mixture ensemble 

1 - 7 exp(- J R 2 g 2 ) 



e(Q) = 0.75 



l-0.57exp(-i? 2 Q 2 



(2) 



with 7 = —2/3. The distribution e(Q) is shown 
in Fig. [2| fitted with the parametrisation given in 
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Figure 2. The fraction e(Q) of the S = 1 con- 
tribution for the (AA, AA) data. The curve rep- 
resents a fit using the parametrisation given in 
Eq. (§). 

Eq. (J2J) with 7 fixed to one. The state S = 1 dom- 
inates for Q > 2 GeV, but it is suppressed for Q < 
2 GeV. The value of the source size estimated 
from e(Q) is R = 0.14 ± 0.09 sta t ± 0.03 sys fm, in 
agreement with the value obtained from the cor- 
relation function. From comparison with results 
from 7T^~ 7T^~ and K^K^ correlation measurements, 
one observes that the source dimension decreases 
with increasing mass of the emitted particle. 

3. Colour reconnection in W-pair decays 

The colour reconnection in e + e~ — > W + W~ 
was studied in a data sample of 174.2 pb _1 
at a centre-of-mass energy of ^/s = 189 GeV. 
Hadronic and semilcptonic W-pair decays were 
selected and the experimental distribution — In x p 
of the charged particles was compared for each 
event type to the MC models KORALW and 
EXCALIBUR without CR and to EXCALIBUR 
with CR. Here x p = p/(y / s/2) is the scaled mo- 
mentum of a particle. KORALW and EXCAL- 
IBUR are used to generate W-pairs; both are 
coupled to JETSET for the hadronization part. 
Three CR models §, denoted as I, II and IIP 
and implemented in JETSET, were compared to 



data. The distribution obtained for the data 
and the MC models are shown in Fig. ||. The 



4. Bose Einstein correlations in W-pair de- 
cays 
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Figure 3. The — \nx p distributions for semilep- 
tonic and hadronic W-pair decays for data and 
MC models with and without CR. 

ratio of the multiplicity in fully-hadronic events 
to twice the multiplicity in semileptonic events 
is also shown. The multiplicities within the ex- 
perimental acceptance for the semileptonic chan- 
nel are N' 
and NjZ* 



ch 



= 17.53 ± 0.19 ± 0.24 for data 
17.41 ± 0.04 ± 0.29 for MC (no 
CR), giving a difference between data and MC 
of 0.12 ±0.42. For the fully hadronic channel, the 
multiplicity is iV™ g = 35.52±0.22±0.43 for data 
and = 34.77±0.04±0.58 for MC (no CR), 

which gives a difference of 0.75 ±0.76. The differ- 
~ is 0.47±0.44± 0.26 for data, 



enceA^ qqq - 



-2N 



Zz^qq 



ch 



-0.05 ±0.09 for MC (no CR) and 0.52 ±0.52 for 
the difference between data and MC. No colour 
reconnection was observed in the data, but at the 
current level of statistical precision both models 
with and without CR are compatible with the ex- 
perimental results. 



The Bose-Einstein correlations in W-pair de- 
cays were studied using data recorded by the 
ALEPH detector at centre-of-mass energies of 
172, 183 and 189 GeV. For the tuning of the MC 
models of BE correlations, Z data recorded at 
91.2 GeV with the same detector configuration 
as for W + W~ events were used. Pairs of unlike- 
sign pions were chosen as reference sample and 
the correlation function was defined as 
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The correlation function was parametrised with 
R*{Q) = K (l + eQ)(l + \c-° 2 Q 2 ) (4) 

where the term 1 + Ae _cr y describes the BE ef- 
fect. The two parameters A and a characterise 
the effective strength of the correlations and the 
source size, respectively. The term 1 + eQ takes 
into account some long range correlations due to 
the charge and the energy-momentum conserva- 
tion, while k is a normalisation factor. 

The BE correlations were first measured in Z 
decays. A MC simulation of the BE effect with 
the JETSET BE 3 model [§ was then tuned on 
these data. The tuned parameters were Ai npu t = 
2.3 and i? input = 0.26 GeV. As W bosons do 
not decay into b quarks, the BE effect was de- 
termined separately in an udsc sample and in a b 
sample. Two b samples of different purities were 
tagged (6) and the parameters Ab and A u( jsc were 
determined. Using these parameters, the bb com- 
ponent was replaced by an udsc component. The 
agreement between the MC model and the udsc 
data is good; residual discrepancies between them 
are corrected bin by bin. 

The prediction of the MC model tuned and cor- 
rected on Z data was then checked on semilep- 
tonic W-pairs decays. They were found to be 
in very good agreement. For the hadronic W- 
decays, two cases were considered in the MC 



model: pions from different W's may exhibit BE 
correlations (denoted as BEB) and only pions 
from the same W exhibit BE correlations (de- 
noted BEI). The comparison between these two 
MC models and data is shown in Fig. |J. The re- 
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Figure 4. The correlation function R*(Q) for fully 
hadronic W-decays, compared to MC models of 
BE correlations. 

suit of the fit with the parametrisation given in 
Eq. (^) is also shown. All four parameters k, e, A 
and a were free in this fit. The values of A and a 
were used to compute an integral of the BE signal 
Ibe = JiT Xe-" 2 Q 2 dQ oc A/a for data and MC. 
A second fit with k, e and a fixed and A free was 
also made to the first four bins only, where the ef- 
fect is expected to be maximum. The differences 
between the data and the MC models for A from 
the one-parameter fit are 



^data _ ^MC BEB 
\data \MCBEI 



0.088 ± 0.026 ±0.020 
0.019 ± 0.026 ±0.016 



while for the Ibe quantity they are 



rdata 
1 BE ' 

rdata 



I 



MC BEB _ 



BE 



rdata rMC BEI 
1 BE ~~ 1 BE 



= - 0.0217 ± 0.0062 ± 0.0048 
0.0040 ± 0.0062 ±0.0036 



The first error is the statistical error, the second 
is the systematic one. A better agreement is ob- 
tained for the JETSET model with BE correla- 
tions present only for pions coming from the same 



W boson. The JETSET model which allows for 
BE correlations between pions from different W's 
is disfavoured by both A and Ibe variables at 2.7a 
level. 

5. Bose Einstein correlations in W-pair de- 
cays at OPAL 

The OPAL collaboration has analysed data 
recorded at centre-of-mass energies of 172, 183 
and 189 GeV. Three mutually exclusive event 
sample were selected: the fully hadronic event 
sample W + W~ — > qqqq, the semileptonic event 
sample W + W _ — > Z^qq and the non-radiative 
(Z°/7)* event sample (Z°/-f)* — > qq. The cor- 
relation function C(Q) was defined according to 
Eq. (||). For each sample, the correlation func- 
tion was written as a combination of contribu- 
tions from the various pure pion classes, includ- 
ing the background. For hadronic event sample 
one has 

C had (Q) =P^T(Q)C^(Q) + 



WW/ 



had (Q)]C& (Q) 



Similar expressions were written for C seml (Q) and 
C non - rad (Q) correlation functions. The prob- 
abilities Ph[]^(Q), etc. were taken from MC 
simulations without BE effect. Each correla- 
tion function C* q ™(Q), C™(Q) and C z * (Q) was 
parametrised by 

C{Q) = N[l + U(Q)XeM-Q 2 P 2 )} (5) 

where /^(Q) is the probability of the pair to be 
a pair of pions. A simultaneous fit was made to 
the experimental data, with the parameters N, A 
and R free for each event class (nine free param- 
eters). All three classes exhibit BE correlations 
with consistent R and A parameters. BE correla- 
tions were then investigated separately for pions 
coming from the same W and from different W's. 
The correlation function for the hadronic event 
sample was written as 

c had (Q) = p£r(Q)c samc (Q) + p h z ; d (Q)c fc z ; cq) 

+ [i - P^lTiQ) P h z ; d (Q)]c diff (Q) 

where C samc (Q), C* diff (Q) and Cf g {Q) are the 
correlation functions for pions from the same 



W, from different W's and from (Z°/7)* — * qq 
events. Similar expressions were written for 
C scmi (Q) and C non - rad (Q). The correlation func- 
tions C samo (Q), C diff (Q) and C Z *(Q) were un- 
folded from the data and are shown in Fig. |[ 
They were then parametrised by Eq. (^|) and si- 
multaneous fits to the experimental distributions 
were performed. Three different cases were con- 
sidered: 1) the same source-size R for all event 
classes; 2) different R parameters for each class, 
and 3) i? diff is related to i? same using the theoret- 
ical prediction (i? diff ) 2 = (R^™) 2 + 4/3 2 7 2 c 2 r 2 . 
The results obtained for the parameter A for the 
third case are 

A same = 0.69 ± 0.12 ± 0.06 

A diff = 0.05 ±0.67 ±0.35 
A z * = 0.43 ±0.06 ±0.08 

At the current level of statistical precision it is 
not possible to determine if correlations between 
pions from different W's exist or not. 
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Figure 5. The correlation function for unfolded 
classes C samc (Q), C diff (Q) and C Z *(Q). 

6. Conclusion 

The ALEPH analyses of Fermi-Dirac correla- 
tions of (AA, AA) pairs in hadronic Z decays and 
of colour reconnection and Bose-Einstein corre- 
lations in W-pairs decays have been presented. 



A depletion of events are observed for the region 
Q < 2 GcV in the (AA, AA) system. In the anal- 
ysis of W-pair decays, no colour reconnection ef- 
fects are observed, but models which predict such 
effects can not be excluded. The Bose-Einstein 
correlations measured in W-pair decays are repro- 
duced by a Monte Carlo model with independent 
fragmentation of the two W's, while a variant of 
the same model with Bose-Einstein correlations 
between decay products of different W's is dis- 
favoured at 2.7cr. The OPAL analysis of Bose- 
Einstein correlations in W-pair decays has also 
been presented. At the current level of statistical 
precision it was not possible to determine if cor- 
relations between pions from different W's exist 
or not. 
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A. Simon, University of Freiburg (Germany) 
The experiments obviously use different methods 
to study BE correlations. It is technically possi- 
ble to agree on one dataset and check the various 
methods on their systematics? 

Answer: Each experiment has optimised the 
analysis on detector characteristics and on the 
aspects of BE correlations which were considered 
important to be studied with the available statis- 
tics. The methods are therefore different; this is 
a normal situation, there is no need to have coor- 
dinated analyses until each experiment obtains a 
set of "final results". A too early coordination and 
convergence of the methods can reduce the qual- 
ity of the analyses. Once the experiments obtains 
"final results", an intercomparison is meaningful, 
followed by cross checks of the methods used by 
the other experiments ( within the limits of statis- 
tics and... manpower). In fact, a LEP group was 
formed recently to understand the differences be- 
tween the results of the four LEP experiments. 



